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Abstract 

We solve a truncated system of coupled Dyson-Schwinger equations for the gluon and ghost 
propagators in SU(A''c) Yang-Mills theories in Faddeev-Popov quantization on a four-torus. 
This compact space-time manifold provides an efficient mean to solve the gluon and ghost 
Dyson-Schwinger equations without any angular approximations. We verify that analytically two 
power-like solutions in the very far infrared seem possible. However, only one of these solutions 
can be matched to a numerical solution for non-vanishing momenta. For a bare ghost-gluon vertex 
this implies that the gluon propagator is only weakly infrared vanishing, Dgi{k^) oc {k'^)'^'^~^, 
K ~ 0.595, and the ghost propagator is infrared singular, Dg/i(/c^) oc (k'^)^'^^^. For non-vanishing 
momenta our solutions are in agreement with the results of recent SU(2) Monte-Carlo lattice 
calculations. The running coupling possesses an infrared fixed point. We obtain a(0) = 8.92/ Nc 
for all gauge groups SU(A'^c)- Above one GeV the running coupling rapidly approaches its 
perturbative form. 
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I. INTRODUCTION 



It is generally accepted that the theory of strong interactions, QCD, should describe the 
observed phenomenon of confinement: colored objects like quarks and gluons occur only in 
hadrons. A possible route for gaining more understanding of this phenomenon is the study 
of the infrared behavior of QCD Green's functions, for a recent review see [|I|. In addition to 
shed light on the fundamental properties of QCD the knowledge of these Green's functions 
provides the basis for a successful description of hadronic physics ^ 0]. Based on the idea 
of infrared slavery older works on this subject assumed a strongly infrared singular gluon 
propagator. Recent studies based either on Dyson-Schwinger equations 0, ||, ||, |^, or 
lattice calculations ^ [1^, |TT| in Landau gauge indicate quite the opposite: an infrared 



finite or even infrared vanishing gluon propagator. These two techniques are complementary 
in the following sense: On the one hand, Monte-Carlo lattice calculations include all non- 
perturbative physics of Yang-Mills theories but cannot make definite statements about the 
very far infrared due to the finite lattice volume. On the other hand, Dyson-Schwinger 
equations allow to extract the leading infrared behavior analytically and the general non- 
perturbative behavior with moderate numerical effort but these equations, consisting of an 
infinite tower of coupled non-linear integral equations, have to be truncated in order to be 
manageable. As we will also see in the course of this article, the propagators of SU(2) and 
SU(3) Landau gauge Yang-Mills theory in Faddeev-Popov quantization coincide for these 
two different approaches reasonably well. Thus we are confident that our results for the 
qualitative features of these propagators are trustable. 

Especially, these recent results on the Landau gauge propagators imply that the Kugo- 



Ojima confinement criterion |T2|, |T^, |I^ is fulfilled (for a short summary on this topic see e.g. 
ref. [0). It is gratifying to note that no truncation to the Dyson-Schwinger equations has 
to be applied to arrive at this conclusion if one assumes that the involved Green's functions 
can be represented in the infrared by asymptotic expansions 0, 0. In Landau gauge, a 
sufficient condition for the Kugo-Ojima confinement criterion is that the nonperturbative 



ghost propagator is more singular than a massless pole in the infrared |I^ : 



D^\p) = -S'^'^ , with G{p')'^ oo. (1) 



This behavior is also correlated to other aspects of Yang-Mills theories. First, the Oehme- 



Zimmermann superconvergence relations |16| can be derived from Ward-Takahashi identities 



assuming the Kugo-Ojima confinement criterion ||T^. These superconvergence relations 
formalize a long known contradiction between asymptotic freedom and the positivity of 
the spectral density for transverse gluons in the covariant gauge. Second, eq. (0) agrees 



with Zwanziger's horizon condition |^, |T8[. This amounts to Gribov's prescription to cut 
off the functional integral at the first Gribov horizon [|T^. Noting that this horizon is a 
convex hypersurface in A-space that surrounds the origin [EO] allows one to conclude that 



the Dyson-Schwinger equations are not changed 113. However, one has to note that this 



treatment of the functional integral is related to the resolution of an ambiguity in the solution 
of these equations |^. Nevertheless, supplementing the Faddeev-Popov quantization with 
this additional constraint might be not sufficient to provide an exact solution of the problem 
because there exist Gribov copies within the first Gribov horizon . 



As already stated, the obtained values for the infrared exponents of the gluon and ghost 
propagators depend on the employed approximation for the Dyson-Schwinger equations. 
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Beyond the necessary truncation of this set of integral equations in numerical calculations 
also some approximations for the angular integrals have been used so far Q . On the other 
hand, employing infrared expansions (without using any angular approximation) which are 
strictly valid only in the limit of vanishing momentum, 0, 0, |^, ^, yield also 

some quite different values for the infrared critical exponent depending on the truncation 
scheme. With respect to these analytical calculations the question arises whether for every 
extracted value of the infrared exponent the corresponding numerical solution exists also 
for finite values of momenta. As will be detailed in this paper a tool to overcome angular 
approximations is the treatment of the Dyson-Schwinger equations on a compact Euclidean 
four-manifold, in our case a four-torus. This allows us to answer the above question: We will 
see that not every analytically extracted infrared exponent can be matched to a numerical 
solution for non- vanishing momenta. 

To make this paper self-contained we will shortly summarize truncation schemes for the 
gluon and ghost Dyson-Schwinger equations in fiat (infinite) Euclidean space-time, which 
have been solved recently, in the first two sections of Sect. II. In the following section we 
introduce a novel truncation scheme. In Sect. Ill we present the Dyson-Schwinger equations 
formulated on the momentum grid which is the dual space to the compact four-torus. In 
Sect. IV we will present solutions done in otherwise exactly the same approximation scheme 
as previous solutions for fiat space-time. The comparison to these previous solutions allows 
to chose a suitable regularization and renormalization procedure. (As in previous work 
we adopt a modified MOM scheme.) In Sect. IV also the solutions without any angular 
approximations will be discussed. The central result of this paper is: only one of two 
solutions allowed by the infrared analysis can be matched to a numerical solution for non- 
vanishing momenta. E.g. for bare vertices this implies that the gluon propagator is only 
weakly infrared vanishing, Dgi{k^) oc (/c^)^''"^, /t = 0.595..., and the ghost propagator is 
infrared singular, Dgh{k'^) oc {k'^)^'^^^. In Sect. V we present our conclusions. Furthermore, 
in Appendix A we present our approximation to impose one-loop scaling of the gluon and 
ghost propagators in the ultraviolet. In Appendix B we discuss the influence of gluon and 
ghost zero modes on the solutions. The numerical methods to solve the gluon and ghost 
Dyson-Schwinger equations in fiat Euclidean space-time will be described in Appendix C. 

II. GLUON AND GHOST DYSON SCHWINGER EQUATIONS IN FLAT EU- 
CLIDEAN SPACE-TIME 

In this section a short summary of previously employed truncation and approximation 
schemes for the coupled gluon and ghost Dyson-Schwinger equations 0, Q will be given first. 
We will also provide the underlying formula for a new truncation scheme. All these schemes 
include all diagrams in the ghost equation and neglect contributions from the two-loop 
diagrams in the gluon equation, see Fig. |l], where the full gluon Dyson-Schwinger equation 
of QCD is represented diagrammatically. In addition, as we will be only concerned with 
pure Yang-Mills theory in this paper the quark loop will be neglected. The tadpole term 
provides in Landau gauge only an (ultraviolet divergent) constant and will drop out during 
renormalization anyhow. Thus, we will effectively study the coupled system of equations as 
depicted in Fig. |^. 

The necessary trunction of the gluon Dyson-Schwinger equation immediately leads to 
a problem: The gluon polarization, which due to gauge symmetry would be transverse to 
the gluon momentum in an exact calculation, acquires spurious longitudinal terms. These 



3 



FIG. 1: Diagrammatic representation of the gluon Dyson-Schwinger equation. The wiggly, dashed 
and sohd hnes represent the propagation of gluons, ghosts and quarks, respectively. A filled blob 
represents a full propagator and a circle indicates a one-particle irreducible vertex. 
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FIG. 2: Diagrammatic representation of the truncated gluon and ghost Dyson-Schwinger equations 
studied in this article. In the gluon Dyson-Schwinger equation terms with four— gluon vertices and 
quarks have been dismissed. 
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terms are in general quadratically ultraviolet divergent and thus highly ambiguous because 
they depend on the momentum routing in the loop integral. In addition, a gauge invariant 
regularization scheme is required to avoid these unphysical longitudinal terms. Such schemes 
are, however, hard to implement in Dyson-Schwinger studies, for the corresponding use of 
dimensional regularization see e.g. refs. [^, An alternative unambiguous procedure is 
to isolate the part free of quadratic ultraviolet divergences by contracting with the projector 

n,^{k) = V - = V - 4^ , (2) 

which is constructed such that lZ^u{k)5^^ = 0, and therefore the ambiguous term propor- 
tional to 6^^ is projected out ||2^. Note that the use of this projector also removes the 
tadpole term. As has become obvious recently (see also ref. for a corresponding dis- 
cussion in a much simpler truncation scheme) the use of the projector interferes with the 
infrared analysis of the coupled gluon-ghost system. For technical reasons we will employ a 
one-parameter family of projectors 

'Pli\k)=5,.-C^, (3) 

which allows us to interpolate continuously from the projector (Q) to the transversal one 
(with = 1). Furthermore, we will also use the general form of the quadratically ultraviolet 
divergent tadpole to remove these unwanted ultraviolet divergencies. 

Here we employ the conventions and notations of ref. As usual for Dyson-Schwinger 
studies all integrals are formulated in Euclidean space-time. Z{k'^) is the gluon renormal- 
ization function defined via the gluon propagator in Landau gauge 

D^M = i-* (s,. - ^) ^ . (4) 

The deviation of Z(k'^) from its tree-level value Z = 1 provides a measure for renormalization 
of the gluon field due to the considered interactions. The ghost renormalization function 
GikP') is defined analogously via the ghost propagator, see eq. (|1]). 

In Landau gauge the ghost-gluon vertex does not attribute an independent ultraviolet 
divergence, i.e. one has Zi = 1, [^. Therefore a truncation based on the tree-level form for 
the ghost-gluon vertex function, G^{q,p) = iq^ is compatible with the desired short distance 
behaviour of the solutions. This will be exploited in the following. 

To proceed we will first consider the truncation schemes of refs. ||^, ^ in the next two 
subsections. Both these truncation schemes employ only the projector lZi^,y (^ in the gluon 
equation. The main difference between these truncation schemes consists in the treatment of 
the three-point functions. Whereas in ref. |^ the form of ghost-gluon and three-gluon vertex 
function has been related to the gluon and ghost renormalization functions using Slavnov- 
Taylor identities (and then the resulting system has been solved self-consistently), in ref. Q 
bare three-point functions have been used. Amazingly, though, both schemes provide results 
with identical qualitative infrared behavior: the gluon propagator vanishes in the infrared, 
the ghost propagator is highly singluar there, and the strong running coupling (which can 
be related to the gluon and ghost renormalization functions using the specific form of the 
ghost-gluon vertex in Landau gauge 0) has an infrared fix point. Because this infrared 
behavior is determined by the interplay between the ghost loop in the gluon equation (the 
gluon loop being subleading in the infrared) and the ghost equation such a scenario is also 
found in the ghost-loop only approximation We will exploit these two recent truncation 
schemes in the course of this article for various tests of our method. 
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A. The dressed vertex truncation including the gluon loop 



In ref. an approximation scheme for the longitudinal parts of ghost-gluon and 3-gluon 
vertex functions has been employed which ensures consistency at the level of one-particle 
Greens functions, i.e. propagators. The detailed form of the vertex functions can be found 
in ref. Q]. Here we provide directly the ghost, 



p 



and the gluon 



1 _ ^.«!^/ A f^^,^._^.,.^(pW)^tf)Gtf) 



...2 2 .2^Z{p^)G{p') ,^,2 2 .2^Z{q^)G{q^)\ G{k^) 

+^ / (^(^')^(P') - G{q')G{p')) - {Qnk)p)G{e)G{q')^^^^ , 

equations, respectively. The functions Ni{x,y;z) = Ni{y,x]z) and N2{x,y;z) are given in 
appendix C of ref. 0. 

In eq. (|^) we have already exploited the identity Z\ = \. This leaves the gluon and 
ghost field renormalization constants Z3 and Z3 as well as the gluon vertex renormalization 
constant Z\ to be determined correspondingly to the employed truncation. Note that these 
constants depend on the ultraviolet cutoff A and the renormalization scale /i. 

In ref. different angular approximations for > k^ and for < k'^ have been em- 
ployed. In the latter case G(p^) = G[(k, — g)^) G(k^) and Zij?') Z{k'^) have been set 
which obviously preserves the limit g^ — of the integrand. With this approximation one 
obtains from eq. @ upon angular integration 



G{k^) " 167r2 4 Wo 




Z{e)G{e)+ [ %Z{q^)G{q' 



^3 - 1 1 z{e)G{k^) + r ^4 1 , (7) 



2 ifc2 q 

where we introduced an 0(4)-invariant momentum cutoff A to account for the logarithmic 
ultraviolet divergence, which will have to be absorbed by the renormalization constants. 
With some further assumptions the angular approximation for the gluon equation (|^) yields 



Z{k^) 'IQtt^ 3 I 7o V2A;4 2 P 

dq^ [Ik 



q2 



167r2 3 Wo A;2 2 A;2 ^ ' ' 3 ^ ' 2 J^2 q^ 
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In the infrared the solutions Z{x) and G{x) behave power-like, 



Z{x) oc G{x) oc x-\ (9) 

In this truncation scheme one obtains k ~ 0.92 0. To solve the coupled system for all 
momenta the power laws, eq. (^, are used to perform the integrals from ?/ = to an 
infrared point y = e analytically, while the remaining part of the integrals is done with the 
help of numerical routines, see e.g. ref. 



The non-perturbative subtraction scheme of ref. ^ implies a strong running coupling 
with infrared fixed point. Starting again from the non-renormalization of the ghost-gluon 
vertex _ _ 

Zi = ZgZlJ^Z^ = 1 , (10) 

one can readily show that the product g'^Z{fi^)G'^{fi'^) is renormalization group invariant. 
Therefore, in absence of any dimensionful parameter, this (dimensionless) product is a func- 
tion of the running coupling g only. Analysing this renormalization group invariant product 
more closely one concludes that it is identical to the running coupling g'^{ii^). As the in- 
frared powers in the product ZG^ cancel exactly the running coupling is a finite constant 
for /i^ = 0. 



B. The bare vertex ghost-loop only truncation 

Substituting the tree-level ghost-gluon vertex for the dressed one and neglecting the gluon 
loop the coupled system of equations (^,^ reads 

' - Z3 - g-NX [ ^ ff~!^'f G{q^)Z{{k - qn (H) 



1 _ ^N^~ r dS k\k^ + + {k - qf) + Aq\k - qf - 2{k - q)^ - 2q'^ 

'ZiJ^) ~ '^^TV(2^ 2k^q^{k-qy 

xGiq')Giik-q)'). (12) 

To simplify notation we introduce the abbreviations x := k'^ , y := q"^ , s := fi"^ and L := h?. 
In the angular approximation of ref. ^ ('ymax approximation') where Z(min(p^, g^)) and 
G(min(p^, q^)) are substituted for Z{{k — g)^) and G{{k — qY) eqs. (|lT],|12D are simplified to 



Imposing as renormalization conditions Z{s) = G{s) = 1 to determine Z^ and Z3 these 
equations may be solved numerically for a given (0(4) invariant) cutoff. As a matter of fact, 
we use subtracted finite equations in our numerical procedure, see Appendix C for details. 



In the infrared, also the solutions of eqs. ([T3|JT^ ) behave power-like, c.f. eq. (|), with 



K ^ 0.77 0. In the same truncation but with no angle approximation employed a solution 
K = 1 has been extracted for the infrared behaviour. However, it will be explained below 
that we could not find a numerical solution for non-vanishing momenta connected to the 
K = 1 infrared behaviour. 
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C. The bare vertex truncation including the gluon loop 



In this subsection we will detail a novel truncation scheme which employs bare vertices. 
Nevertheless it will be constructed such that it reproduces the correct perturbative limit for 
large momenta. To analyse the gluon loop we will use the class of projectors V^J (|). A 
smooth interpolation between the Brown-Pennington projector {( = 4) and the transverse 
one {( = 1) will be helpful in the analysis of unphysical quadratic ultraviolet divergencies 
for 7^ 4 in the gluon equation. Their careful removal is essential for a numerical solution 
retaining the infrared behaviour of the solutions. 

The coupled equations for the ghost and gluon dressing functions using bare vertices read 
as follows (The ghost equation is, of course, identical to eq. (|ll|), it is only repeated for a 
coherent represent at ion . ) : 

! #^^^G(,)Z(.) , (15) 



G{x) ^ i (27r)4 xy 

1 7 , 2^0 f d^q M{x,y,z) x , ^ f A Q{x,y,z) 

- Zs + g — ——- G{y)G{z) + Zig — / — — Z{y)Z{z) 



Z{x) ^37 (27r)4 xy ^ 3 J (27r)4 xy 



The kernels ordered with respect to powers of z := p = {k — q) have the form: 



(16) 



r.( X I f {x-yf\ 1 fx + y\ 1 



M{x,y,z) = -{——x + ---—] +- + ----- (18) 



^^4 2 4 X J 2 2x Ax 

1 f^x^ 2 19-C 5-C 2 Cy^ 

Q{x,y,z) = + x — xy + ——y + - — 

z'^ \8 y 8 4 8x 

1 fx' 15+c 17-c 

+ -{ -. — X — y + C— 

z \ y 4 4 X 

19-C a: ^ 17-C , 9C ^ 
8 y 4 Ax 

+ (19) 

\x Ay J Sxy 

It is straigtforward to verify that for C = 4 the kernel M{x, y, z) is identical to the kernel in 
eq. (|12D. 

The quadratic ultraviolet divergencies of the integrals are most easily discovered approx- 
imating the angular integrals as done in the previous subsections (note that we will use 
these approximations only for an analysis of the ultraviolet behaviour) and introducing an 
Euclidean sharp cutoff. Displaying only the ultraviolet divergent integrals 

r ? ((^ - -^<'^) - ^ 

one sees that for C = 4 the terms independent of the integration momentum y vanish. Thus 
these integrals are then only logarithmically ultraviolet divergent as could be expected on 
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the basis of the results summarized in the previous subsections. Of course, the quadratic 
ultraviolet divergencies are artefacts of the employed truncation. Due to gauge invariance 
they would cancel against the tadpole and similar divergencies in the two-loop terms. The 
calculation of the latter being beyond the scope of this paper we will simply subtract these 
divergent terms. This cannot be done straightforwardly at the level of integrands: Such 
a procedure would disturb the infrared properties of the Dyson-Schwinger equations. As 
we anticipate from previous studies and analytic work [^] that the ghost loop is the 
leading contribution in the infrared the natural place to subtract the quadratically ultraviolet 
divergent constant is the gluon loop. We do this by employing the substitution 



Q{x,y,z) Q'{x,y,z) 



Q{x,y,z) + -{A-C) 



(21) 



in eq. (|TB|). At the first sight, due to the presence of the prefactor Zi in the gluon loop, this 
seems not to be sufficient to remove also the quadratic ultraviolet divergence of this ghost 
loop. However, note that in the next step we will enforce consistency of the logarithmic 
divergencies which entails then cancelation of the quadratic divergencies if Q' is employed 
in the gluon loop, see below. 

To achieve the correct one-loop scaling in the ultraviolet we will adopt a similar treatment 
to the one of ref. 0]. Please note that within the presented class of truncation schemes 
it is impossible to fulfill both, correct one-loop scaling and the Slavnov-Taylor identity 
Zi = Z3/Z3. For large Euclidean momenta and to one loop the behavior of the propagator 
functions can be described as 



Z{x) = Z{s) 
G{x) = G{s) 



log ( ^ ) + 1 
ulogi^] +1 



(22) 
(23) 



Z{s) and G{s) denote the value of the dressing functions at some renormalization point 
s := yU^, 7 and 6 are the respective anomalous dimensions. To one loop one has 6 = —9/44 



and 7 = — 1 — 25 for arbitrary number of colors Nc and no quarks, Nf 
u = llNca{s)/127i. 







Furthermore, 



Employing these expressions in the ghost equation ([TSD and approximating the angular 
integrals as done previously one obtains to the order of the leading logarithms 



9g'Z{s)G{s) 



X 







7+5+1 






ulog (^j^ +1 




^log +1 











-y+S+l 



(24) 



where the abbreviation L := A^ has been used again. The renormalization constant Z^{L, s) 
cancels the cutoff dependence of the first term in the bracket. Thus, the power and the 
prefactor of the second term have to match with the left hand side of this equation. This 
leads to two conditions: 



7 + 25 + 1 
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9 ^,_^^2/ 



2t,;(7 + 1) levr^ 



;Zis)G' 



(25) 
(26) 
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Eq. (p5D is of course nothing else but consistency of the ghost equation with one-loop scaling. 

Plugging the ultraviolet behavior (^2|) and (^) into the gluon equation (0) and keeping 
the leading order logarithms leads to 



ijJ log 






(^) 


+ 1 









327r^ct;7 
culogi J 1 +1 



lg'Z{s)G{s) 



dyZx 



log ( - 

s 



wlog(^) +1 



-7 



7+<5 



(27) 



Here we see that due to the employed truncation we have to give up either one-loop scaling or 
the Slavnov-Taylor identity Z\ = Z3/Z3. Note, however, that Zs{L, s) has the correct cutoff 
dependence to make this equation finite with finite Zi{L,s). Nevertheless, in order to get 

the correct leading log contribution from the gluon loop the integrand in the last line should 

-7-1 



be proportional to 



uj log 



+ 1 



. This not only contradicts the identity Zi = Z^/Z^ 

but also requires the renormalization constant Zi{L, s) to acquire a momentum dependence 
if correct one-loop scaling is to be enforced. Following case c) described in chapter 6 of ref. 
|3ll we demand therefore 



Zi{L,s) Zi{x,y,z;s,L) 
Generalizing the approach of ref. 

Zi{x,y,z;s,L) 



log ( - 



-5-27-1 



(28) 



we note that the following two-parameter ansatz 

^(^)(l-a/5-2a)^(^)(l-6/5-26) 



Z{y)(^+^) 



(29) 



with arbitrary a and b satisfies the required proportionality. This can be straightforwardly 
verified with the help of expression (|25|) obtained from the ghost equation. In addition 
ansatz (^91) ensures the cancelation of quadratic ultraviolet divergencies as discussed above 
also if C 7^ 4. 

Carrying out the remaining integral in (|27|) and matching powers one obtains 



9' 



2'yuj IGn' 



-Ms)G\s) 



:Z{s 



7 + 25+1 

bQ^^Y-a/5-b/S-2a-2b 



0, 

1 . 



(30) 
(31) 



First of all, we reproduce correct one-loop scaling also from the gluon equation as can be 
seen from the equivalence of eqs. (0) and (^5]). Second, with the perturbative renormaliza- 
tion conditions Z{s) = G{s) = 1 which are possible and appropriate at a large Euchdean 
renormalization point /i one obtains together with eq. (^) the correct anomalous dimensions 
7 = -13/22 and 6 = -9/44. 

A reasonable choice of parameters is of course one which keeps Zi as weakly dependent 



as possible on the momenta y and z, c.f. Fig. |10] in Appendix A. Note that the choice 
a = b = corresponds to the truncation scheme of p8| whereas a = 36,b = together 
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with the appropriate vertex dressings reproduces case c) of ref. 0. The infrared behavior 
of the gluon loop in the gluon equation depends strongly on a and h. E.g. for 6 = one can 
distinguish three cases: For a < the gluon loop is subleading in the infrared, for a = as 



in ref. |^ the gluon loop produces the same power as the ghost loop, for a > the gluon 
loop becomes the leading term in the infrared. In the latter case we did not find a solution 
to the coupled gluon-ghost system. In Appendix A we will demonstrate that a = = 35 
minimizes the momentum dependence of Zi. Thus we will use these values except stated 
otherwise explicitely. 

The infrared behavior of the propagator functions in this truncation scheme and for the 
transverse projector {Q = 1) has been determined very recently i, 0. It can be shown 
straightforwardly that the relation between ghost and gluon infrared behavior is again as in 
eq. (P). This is expected also on general grounds f^. Due to this one can immediately use 
the ansatze 

Z{x) Ax^^ , G{x) Bx-'' for x^O . (32) 

As is explained in detail in ref. [0 the renormalization constants and Z3 can be dropped 
for very small momenta x: They are either subleading in the infrared (gluon equation) 
or have to be zero when the renormalization takes place at /i = (ghost equation). The 
remaining infrared integrals can be evaluated using the formula |^ 



. 2+a+6r(2 + a)r(2 + 6)r(-a-6-2) 



r(-a)r(-6)r(4 + a + 6) ' ^ ^ 



where again x = k'^, y = q'^ and z = {k — qf. From the resulting two conditions 

\AB^ , (34) 



1 (2 + «:)(! + «) _ T\2 - k)T{2k) g^N, ^ . 



18 (3-2fi:) r(4-2K)r2(l + K) 487r2 

4fi;-2 T^{2-k)T{2k) g^N, 



4C«;-4k + 6-3C T{A-2k)T^{1 + k) A8n^ 



AB^ , (35) 



one obtains the determining relation for k by equating both left hand sides, see Fig. ^. For 
the Brown-Pennington projector, i.e. ( = 4, one then finds the known solution k = 1 [Q. 
However, as can be seen immediately the left hand side of the second equation possesses a 
zero for k, = 1/2 which is canceled by a pole only for ( = 4. Lowering ( only slightly a 
further solution with k slightly larger than 0.5 exists. For the transverse projector, i.e. ( = 1, 
this latter solution becomes k = 0.59... in accordance with refs. |^. Also the solution 
K, = 1 changes continuously when lowering (. The corresponding k, are then all larger than 
1 and contradict the masslessness condition, see chapter 5 of ref. ^ for a discussion of this 
condition. The main result of this paper is that the infrared behavior k, 0.5 matches to a 
numerical solution whereas no numerical solution could be found with the infrared behavior 

K = l. 

The renormalization group analysis for the running coupling given in ref. certainly 
applies here as well. The renormalization group invariant expression for the running coupling 
is therefore given by 

a{x) = —Z{x)G^{x) . (36) 
An 



As can be seen directly from equations ( p4| , p5D , in the here presented truncations the product 
Ncg'^AB'^ is constant for a given k. With a{0) = g'^AB'^/An one concludes immediately that 
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FIG. 3: Here the graphical solution to equations (|3^ ) and (pq) is shown. The thick line represents 
the left hand side of equation (34), whereas the other four curves depict the left hand side of 
equation ( |35| ) for different values of the parameter The ellipse marks the bulk of solutions 
between k = 0.5 and k = 0.6 for different whereas the circles show the movement of the solution 
K = 1.3 for a transverse projector to k = 1 for the Brown-Pennington case, C = 4. 



a{x) is proportional to A^""^. Furthermore, the ghost and gluon dressing functions Z{x) and 
G{x) are independent of the number of colors: enters the Dyson-Schwinger equations 
only in the combination g'^Nc at our level of truncation. From the solution k, = 0.5953 of 
the infrared analysis with the transverse projector C = 1 o^ig determines the infrared fixed 
point of the running coupling to be a(0) = 2.972 for Nc = 3. 



III. GLUON AND GHOST DYSON SCHWINGER EQUATIONS ON A FOUR- 
TORUS 

There are three central aims connected to the investigation of the Dyson-Schwinger equa- 
tions on a four-torus. The first one is purely technical: This allows to study finite volume 
effects also in Dyson-Schwinger approach. Monte-Carlo simulations on a lattice necessar- 
ily have to be done in a finite volume. Therefore in the latter kind of approach infrared 
properties are only accessible by extrapolations to an infinite volume where the avalaible 
data are gained on several different volumes which, due to limitations in computer time, do 
only differ at best by one order of magnitude. We will see in the present Dyson-Schwinger 
approach that available volumes cover several orders of magnitude. And more importantly, 
in several truncations and approximations one can compare to the results obtained in an 
infinite volume. 

The second issue is the solution of the (truncated) Dyson-Schwinger equations without 
any angular approximations. We will detail below why and how this can be achieved on 
the momentum grid dual to the four-torus. This will also lead us to the main result of this 
paper discussed in the following section. 
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The third aim relates to the possibility of topological obstructions on a compact manifold. 
It is well known e.g. that a four-torus allows for a non-vanishing Pontryagin index. We hope 
that we will be able to describe the solutions of Dyson-Schwinger equations with twisted 
boundary conditions ||2^ in a subsequent publication. And choosing an asymmetric four- 
torus might allow the introduction of a non-vanishing temperature in a relatively simple 
way. 



A. Finite volume effects 

From a technical point of view using as underlying manifold a four-torus or choosing 
periodic boundary conditions on a hypercube is identical. Note that the definition of the 
Faddeev-Popov operator necessitates periodic boundary conditions for ghosts instead of 
their Grassmannian nature, see e.g. ref. [^. With L being the length in every direction 
of the hypercube the four-dimensional momentum integral has to be substituted by a sum 
over four indices, 

On the other hand, the quantities of interest, the gluon and ghost renormalization function, 
and G{k'^), resp., do only depend on the 0(4) invariant squared momenta as long 
as all directions are treated on an equal footing on the torus. This suggests to relabel the 
points on the momentum grid not according to a Cartesian but a hypersherical coordinate 
system, 

h E =^E. (38) 

where the index j numbers the hyperspheres = const. The index /, which numbers the 
grid points on each hypersphere respectively, will be dropped in the following. 

In the integrals to be discretized there appear three momenta, the external momentum, 
labeled k, the loop momentum q and for the second propagator in the loop p = k — q. We 
will use the following notation: 

X := k"^ with Xi E hypersphere i , 
y := (f with Uj G hypersphere j , 
z := = {k — qY with G hypersphere n . (39) 

On the hypercubic momentum grid dual to the four-torus the momentum p = k — q is located 
on the grid for every pair of grid momenta k and q as can be seen from elementary vector 
operations (or from the analogue of cubic lattices in solid state physics). 

Note that we have introduced an 0(4) invariant cutoff A into eqs. (0JD. As we will 
see below in more detail a corresponding regularization of the sums over grid momenta is 
required. In a first step we set up a momentum grid with ji, j2, js, j4 = —N, . . . ,0, . . . N. The 
hypercube property of this lattice seem to suggest in the first place to cut off the sums such 
that over complete hypercubes is summed. Such a method, however, breaks 0(4) invariance 
and introduces sizeable numerical errors. (Note that in lattice Monte-Carlo simulations the 
analysis of the resulting data for 0(4) invariance necessitates special kinds of cuts through 
the lattice [^.) As can be seen from Fig. |^ an 0(4) invariant cutoff of the sums necessitates 



13 



FIG. 4: Scetch of the momentum grid dual to the four-torus and the summation over complete 
hyperspheres indicated by fully drawn circles. The hyperspheres depicted by dashed lines are not 
complete due to the numerical ultraviolet cutoff. 



to neglect the "edges" : The sum extends only over the fully drawn hypersheres, and we omit 
the summation over the dashed ones. 

To be able to compare to already known flat space solutions we will first solve the torus 
analogue of the equations with angle approximations, i.e. eqs. (^H) or eqs. (0,0) • To this 



end we start from eqs. or eqs. (|TT],|T2|), respectively. After formulating the corresponding 
four-dimensional integrals as sums over lattice momenta the functional dependence of the 
propagator functions is approximated accordingly to the rules described in the subsections 
II. A and II. B, respectively. We give the explicit expressions for the bare vertex ghost- 
loop only truncation. The corresponding ones for dressed vertex truncation can be derived 
analogously in a straightforward manner, however, as these are quite lengthy we do not give 
their explicit form. 

The Dyson-Schwinger equations in bare vertex ghost-loop only truncation with angle 
approximation read on the torus 



Zsis, L) - g'N' Yl " G(y,)^(max(a:., y,)) , (40) 



G{xi) ^ Xiyj 



3 



y, . Z,{s, L) + ^ ^ " G(y,)G(max(x., y,)) . (41) 

Z[xi) 3 L"^ ^ Xiyj 

As shown in the previous subsection the momentum p = k — q which fulfills Zn = {k — qY 
lies on the momentum grid, and furthermore, p = k — q is determined uniquely from k and 
q for which Xi = /c^ and yj = q^ by using elementary vector operations. Note, however, that 
might be larger than the ultraviolet cutoff even if \/k^ and are not. The actual 
treatment of the sums in eqs. (1^,^) can nevertheless be done straightforwardly as the exact 
algebraic expressions for the kernels are known. 

The main difference between the equations on the torus and the flat space equations is 
the effective treatment in the infrared. The finite volume in coordinate space leads to a finite 
value of squared momentum for the first hypersphere j = 1. Thus one has not to worry 
about possible infrared singularities. On the other hand, there exists zero modes on the four- 
torus. In all the calculations presented here they are neglected, an estimate of their possible 



14 



contribution is given in Appendix B. Note furthermore that for the numerical solution as 
described in Refs. ^ the disadvantage of dealing with infrared singularities has been 
turned into a benefit: The infrared behavior of the propagator functions has been determined 
via an asymptotice series calculating the corresponding coefficients and powers analytically. 
Matching these series at some finite momentum to the functions as obtained from standard 
numerical techniques finally enabled to numerically solve the coupled integral equations. 
Such a complicated process is not necessary employing a torus. However, one anticipates 
already at this level some deviations in the infrared between the solutions obtained in these 
different ways. In the next section the corresponding numerical results will be dicussed. 
They demonstrate that using a torus as infrared cutoff works surprisingly well. 

B. Dyson Schwinger equations with angle integrals on the torus 

The intricacy of the infrared analysis has prevented so far a solution of the coupled 
gluon-ghost system without angular approximations. Using a torus as an infrared regulator 
opens up the possiblity to solve these equations directly without any recourse to an infrared 
asymptotic expansion. As recently arguments have been provided that the use of a bare 
ghost-gluon vertex is not inferior to employing a dressed one |P, |^ we will solve the system 
in the truncation described in sect. II. C. 

With the above discussed replacement of J — > -p the equations ( 
the torus: 

Us, L) - /AT.ij Yl ^^^^^^G(y,)Z(^„) , (42) 
Zsis, L) + E (43) 

+9"^TI E ^^^^^^^G^-^/^-^'^(y,)G^-^/^-^^(.„)Z-(y,)^-^(^n) . 

As already stated ^/i might be larger than the ultraviolet cutoff even if y/x and ^ are 
not. Again the kernels can be nevertheless evaluated straightforwardly, however, if z with 
L < z < ALis the argument of a propagator function we set z = L, i.e. we approximate Z{z) 
or G{z) then by Z{L) or G{L), respectively. Another more elaborate treatment consists of 
matching the corresponding perturbative ultraviolet tail to the function under consideration. 
This has been applied in some cases to test the viability of the method. 

VI. NUMERICAL RESULTS 

A. Comparing torus solutions to previous results in flat space-time 

The solutions of the Dyson-Schwinger equations on the torus can be compared in two 
different ways to the ones obtained in fiat space-time. First, one can use the solutions of 
the "continuum" equations for a certain cutoff A and certain renormalization scale /i to 
provide the values for Z3(yU^,A^) and Z3(/i^,A^) as input for the equations on the torus 
which are then solved, c.f. the discussion in Appendix C where the numerical methods 
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FIG. 5: Shown are the ghost dressing function, the gluon dressing function and the running cou- 
phng in the bare-vertex ghost-only truncation for different momentum grid spacings corresponding 
to different finite volumes of the torus. The fully drawn lines labeled "continuum" represent the 
respective results for flat Euclidean space-time, i.e. continuous momenta. 



employed to obtain the solutions of the "continuum" equations are summarized. Second, 
one can subtract both equations at the squared momenta sq and sz and therefore trade the 
two renormalization constants for the values of the dressing functions at these momenta, 
namely Z{sz) and G{sg)- For the sake of comparison one can read of these values from the 
continuum solution to be used in the torus equations. If sz and sq are taken to be sufficently 
far in the ultraviolet region of momentum, where finite volume effects play a minor role, the 
two procedures lead to the same results. We have verified that this is indeed the case within 
the limits of numerical accuracy. 

Our results for the ghost dressing function, the gluon dressing function and the running 
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FIG. 6: The same as fig. ^ for the dressed vertex truncation. 



coupling in the bare- vertex ghost-only and the dressed vertex truncation can be seen in Figs. 
^ and ^, respectively. In both truncation schemes we solved for three different momentum 
spacings corresponding to different volumes in coordinate space. To keep the cutoff identical 
for all the spacings within each truncation scheme we have chosen three different grid sizes 
respectively. For the bare vertex truncation they are = 17^, 31^, 51^ and for the dressed 
vertex truncation they are = 21^, 41^, 61^. 

A physical momentum scale cannot be determined in truncation schemes which do not 
provide the correct perturbative running of the coupling in the ultraviolet. With the missing 
gluon loop this is the case in the bare vertex only truncation scheme. We therefore have 
to stick to an internal momentum scale without physical units in this case. The situation 
is different, however, in the dressed vertex truncation scheme. Here we have fixed the 
momentum scale by calculating the running coupling for the colour group SU(3) and using 
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the experimental value a{x) = 0.118 at x = M| = (91.187GeV)^ to fix a physical scale. 

In the bare vertex truncation scheme we have chosen Z{fi^) = G{fi^) = 1 as renormal- 
ization condition where we have taken /i^ = A^, the (squared) ultraviolet cutoff. Of course, 
this choice is by no means special and one is completely free to choose the renormalization 
point wherever one likes. For the numerical calculation in fiat space we have chosen two 
different subtraction points for the ghost and the gluon equations as described in Appendix 
C. For good convergence of the iteration process the ghost equation is most conveniently 
subtracted at zero momentum, whereas the gluon equation can be subtracted at any value 
of squared momentum in the region where the equations are solved numerically. In our 
calculation we chose the cutoff A^ = 0.2 in internal units as subtraction point for the gluon 
equation. This allows us to use the renormalization condition 2'(A^) = 1 directly as input 
in the calculation. The second input is provided by the coefficient A of the leading order 
infrared expansion of the gluon dressing function, ZjFt{x) = Ax^'^. The condition G{h?) = 1 
then leads to A = 357.33. The value of the coupling at the renormalization point, Q;(/i^), 
is taken to be 0.97. Again one is completely free to chose this number arbitrarily up to the 
maximum value of the running coupling which is reached in the very infrared. For the torus 
calculations the unsubtracted equations have been used. The same coupling = 0.97 

has been taken and instead of Z{A'^) and A the values Z^{^'^ = 0.2, A^ = 0.2) = 0.9591 
and Z^{^'^ = 0.2, A^ = 0.2) = 1.1034 for the renormalization constants determined from the 
"continuum" solution have been employed. 

As has already been mentioned, in the dressed vertex truncation scheme we determine a 
"physical" momentum scale according to experimental results. In the numerical treatment 
again the ghost equation has been subtracted at zero momentum whereas the gluon equation 
is subtracted at finite momentum. We solved both equations similar to the method described 
in Ref. especially we introduced also the auxiliary functions F{x) and R{x) as defined in 
Ref. 0. As input values serve the infrared expansion of R{x), R{x) = x'^ + . . ., and the value 
R{s) = 0.8 at the gluon subtraction point s = 1.048GeV^. For the calculations on the torus 
we use the unsubtracted equations with the values Z^{ii^ = M|,A^ = 1.255GeV^) = 1.266 
and Z'i^ij?' = Mf,A^ = 1.255GeV^) = 0.966 for the renormalization constants which have 
been determined from the "continuum" solution. 

Our numerical results, depicted in Figs. |^ and ^, show similar properties for both trunca- 
tion schemes. Compared with the respective "continuum" solutions the ones obtained on a 
torus show deviations for the first few lattice points in the infrared. For higher momenta all 
functions on the torus approach the continuum ones. Deviations are there only visible for 
the curves with the largest spacings. The biggest effect can be seen for the running coupling 
a, which is the "observable" of the system. As a is proportional to the product Z{x)G'^{x) 
the deviations of the dressing functions from the "continuum" curve amplify in the infrared 
in a somewhat erratic way, so that the points in the very infrared cannot be connected by 
a smooth line. Comparing larger and smaller spacings of momentum grids one clearly sees 
that the effect is always one of the first spheres on the respective lattices and therefore moves 
to the infrared for smaller spacings. 

The most important properties of the solutions in flat space can still be found in the torus 
solutions despite some deviations in the infrared. Going from larger to smaller spacings a 
powerlike behaviour of the dressing functions in the infrared with the correct exponents can 
still be infered. For the truncation scheme with dressed vertices the gluon dressing function 
on the torus has the same shape and the same hight of the bump in the bending region of 
the curve. As all curves with sufficently small momentum spacing follow the correct power 
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behaviour within numerical errors we conclude that one surely can adress the question of 
the value of the power k in the infrared on a manifold with finite volume as has been used 
here or as one uses in lattice calculations. This is the central result of the present section: 
Employing a torus as infrared regularization is possible. In the following section we will use 
this to go beyond the angular approximation. 

B. Numerical solutions without angular approximations 

In this section we discuss our numerical results for the bare vertex truncation scheme 
with various projectors in the gluon equation. We show results for the gluon and ghost 
dressing function and the running coupling calculated on a torus without any angular ap- 
proximations. Note that up to now in the gluon-ghost Dyson-Schwinger equations solutions 
for non-vanishing momenta have been obtained only employing some sort of angular ap- 
proximation. Going beyond these approximations we will use our new truncation scheme 
introduced in sect. 2.3. The reason for this is, that the dressed vertex truncation scheme is 
quite complicated, and in its derivation one has dismissed terms which seem to be crucial for 
the infrared behavior beyond the angle approximation. As the survey of gluon and ghosts 
infrared behaviors related to different ansatze for vertex functions [0 does not provide any 
evidence for crucial differences between bare and dressed vertex functions the use of the 
latter seems to be an unnecessary complication. In the case of the bare-vertex ghost-loop 
only truncation with Brown-Pennington projector no solutions for finite momenta without 
angular approximation have been reported yet. In fact the results presented in this section 
suggest that these solutions might not exist. 

This time the values for the renormalization constants and cannot be taken from 
a "continuum" solution. Thus we use the subtracted equations on the torus in the course of 
the numerical solution by iteration. Both equations are subtracted at the renormalization 
point iJ? = 1.9GeV^. Similar to the case of the dressed vertex truncation from the last 
section the 'physical' units are gained from experimental input: a{x) = 0.118 at x = M| = 
(91.187GeV)^. The input values for the dressing functions at the renormalization point are 
Z{iJ?) = 0.83. Requiring G{jj,'^) = 1/ ^/zJjjFj then fixes the overall scale for G{x). For the 
value of the coupling at the renormalization point we chose again a{^^) = 0.97 similar to 
the two truncation schemes in the last section. For the presented solution on three different 
volumes lattice sizes of = 13"^, 43^, 71^ have been used. 

As mentioned in the last section we did check cutoff effects by extrapolating the propaga- 
tor functions at z > A^ with a logarithmic tail with the correct anomalous dimensions. The 
results as compared to the ones obtained by simply setting Z{z) = Z{h?) and G{z) = G'(A^) 
for all z > A^ did change by less than 10"'^. 

Our results for the dressing functions obtained with a transverse projector as shown in 
Fig. 1^ are in agreement with the expected power behavior. The infrared critical exponent as 
calculated in Refs. [§,0]) ^ = 0.5953, thus has been verified: A corresponding (numerical) 
solution for non-vanishing momenta exists. The gluon dressing function is remarkably stable 
against changes of the volume and approaches more and more the expected power solution 
for small momenta. For the ghost dressing function one observes again deviations of the 
first points in the infrared: An extraction of the correct infrared critical exponent from the 
numerical solution for the ghost function is hardly possible. Only some points come close 
to the analytical value of the continuum before the curve starts bending down again in the 
very infrared. For the extracted value of the running coupling in the infrared this leads to 
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FIG. 7: Shown are the ghost dressing function, the ghion dressing function and the running 
couphng in the new truncation scheme without angular approximation for different volumes using 
a transverse projector. 



a distinct mismatch to what one is to expect on the basis of analytical results. 

At first sight the fact that the power solution for the ghost dressing function could be 
not reproduced numerically to a reasonable precision may seem disappointing. Nevertheless 
these numerical results themselves show that the ghost dressing function is highly infrared 
singular. This reflects the long-range correlation of ghosts in Landau gauge. Therefore one 
should expect the ghost dressing function to be the one affected most by a finite volume. 
On the contrary the gluon dressing function vanishes in the infrared and consequently it is 
much less affected by a finite volume. We expect the ghost dressing function together with 
the running coupling to approach more and more the correct power solution in the infrared 
as lattice spacings are decreased and lattice sizes are increased. 
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a(x) 

full angles, ghost and gluon loop, different projectors, SU(3) 




X = [GeV^] 

FIG. 8: The same as fig. |^ for different projectors. 



Furthermore, we add a remark on the transversahty of the gluon propagator as obtained 
in the bare vertex truncation scheme. Of course, the full gluon propagator calculated from 
the complete gluon equation with fully dressed vertices in Landau gauge is transversal due 
to Slavnov-Taylor identities. Thus such a hypothetical solution would be independent of the 
form of the projector, i.e. in our notation with the projector {S^^ — (p^p,^)/p'^ independent of 
the parameter (. In practice, using bare vertices as in our truncation scheme this is certainly 
not the case. Our numerical results for different values of the parameter ( can be seen in 
Fig. p. Although our solutions show the expected dependence on the form of the projector 
this dependence is not too drastic and in general the behaviour of these different solutions 
is very similar. 

For the gluon dressing function one observes that the more ( grows the greater is the de- 
viation from the pure power behaviour. The points in the very infrared cannot be connected 
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FIG. 9: Results on the torus compared to recent lattice results As the torus points are very 
close to each other on a linear momentum scale we did not resolve the torus curves into single 
points. 



by a smooth line any more. Correspondingly this happens for the running coupling. Based 
on the infrared analysis one might anticipate that k should approach the value k = 0.5 more 
and more as C grows until there is a jump to the solution from k = O.S"*" to k = 1 as the 
Brown-Pennington limit C = 4 is reached. We do not observe such a qualitative jump in 
our solutions on a torus. The solution shown for C = 4 is approached smoothly when C, 
approaches this limit. This clearly indicates that the solution k = 1 might not exist at all if 
one removes the torus as a regulator. 

Finally we compare our results to recent SU(2) lattice calculations [^. As has already 
been stated above the ghost and gluon dressing functions from Dyson- Schwinger equations 
are independent from the numbers of colors at least to our level of truncation. The only 
caveat in comparing our results with the lattice ones is the adjustment of the momentum 
scale which is certainly different from the case of SU(3) above. To get an appropriate 
momentum scale we therefore used the lattice result asu{2){'^) = 0.68 at a; = lOGeV^ as 
input. The two graphs in figure (^) show that the main qualitative features, the vanishing 
of the gluon and the divergence of the ghost dressing functions in the infrared, are common 
properties of both, the lattice solutions and the one from Dyson-Schwinger equations. Even 
the power k = 0.59... of the gluon dressing function on the torus is very close to the one 
that can be extracted from the lattice fit to be k ^ 0.5. The main difference of both 
approaches lies in the medium energy region around one GeV, where the Dyson-Schwinger 
solutions suffer from the missing two loop contributions that are certainly present in lattice 
Monte-Carlo simulations. 



V. CONCLUSIONS AND OUTLOOK 



In this paper we have presented numerical solutions of truncated systems of Dyson- 
Schwinger equations for the gluon and ghost propagators in Landau gauge SU(N) Yang-Mills 
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theories. We have employed a four-torus, i.e. a compact space-time manifold, as an infrared 
regulator. This enabled us to overcome angular approximations used so far in previous 
studies The basis for our numerical calculations have been provided by an analytic 

determination |^ of the exponents governing the infrared power-like behavior of the gluon 
and ghost propagators. Typically the infrared analysis provide two possible values each for 
the gluon and the ghost, respectively. The central result of this paper is: Only one of these 
two possible infrared behaviors in a given truncation scheme can be matched to a numerical 
solution for non- vanishing momenta. 

In a truncation scheme with a bare ghost-gluon vertex and a transverse projection of the 
(truncated) gluon equation this implies that the gluon propagator is only weakly infrared 
vanishing, Dgi{k'^) oc (fc^)^''"^, k = 0.59..., and the ghost propagator is highly infrared 
singular, Dgh{k'^) oc (k^)^'^^^. The running coupling possesses an infrared fixed point whose 
value is given by a(0) ~ 2.97 (or, for a general number A^^ of colors, a(0) ~ 8.92/Nc). 

Our results compare very favorably with results of recent lattice calculations performed for 
two colors. Due to the finite lattice volume the lattice results cannot, of course, be extended 
into the far infrared. In this respect our results are complementary to the lattice ones: We 
do obtain the infrared behavior analytically. On the other hand, lattice calculations include, 
at least in principle, all non-perturbative effects whereas we had to rely on truncations. E.g. 
the deviations for the gluon renormalization functions at intermediate momenta depicted in 
Fig. ^ might be due to the neglect of the four-gluon vertex function in our calculations. 

As an outlook we would like to mention that employing a four-torus as underlying mani- 
fold might serve for a number of interesting studies of Dyson-Schwinger equations. Amongst 
these, the use of twisted boundary conditions on the torus might enable one to shed some 
light on the importance of topologically non-trivial gauge field configurations for the infrared 
behavior of QCD Green's functions. We would like to note that recent lattice calculations 



11 1 indicate a relation between the existence of center vortices in the maximal center gauge. 



the area law of the Wilson loop and the infrared behavior of the gluon propagator in Landau 
gauge. 

Furthermore, we want to point out that choosing an asymmetric four-torus might provide 
an efficient mean to extend these calculations to non-vanishing temperatures. Here the main 
qualitative question arises about the fate of the Kugo-Ojima confinement criterion at the 
deconfinement transition. 

Finally, we want to remark that QCD Green's functions are an important input in many 
calculations in hadron physics [l|, The next necessary step towards such phenemenological 
applications is a study of the quark propagator. Such an investigation might hopefully also 
provide some insight into the mechanism of quark confinement in covariant gauges. 
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Z^(y) for different parameter a 




FIG. 10: The y-dependence of the function Zi(y,z) for different values of the parameter a. Only 
the choice a = 35 leads to momentum independence in the infrared. Due to the symmetry of the 
ansatz for Zi{y, z) the z-dependence is the same for b = a. 

APPENDIX A: ONE-LOOP SCALING 

In the framework of the truncation scheme presented in sect. IIC we have shown that the 
approximation 

)(l-a/J-2«)Q( )(l-t/i-2l,) 

Z.(..y..;..L)= y^^^^.^„ (44) 

for the gluon vertex renormalization constant yields the correct one loop scaling of the gluon 
loop in the gluon Dyson- Schwinger equation. This is true for any values a and b. Of course, 
in a full treatment of the coupled ghost gluon system Zi{s,L) would be independent of 
momentum. Therefore a choice of a and b which keeps Zi as weakly varying as possible 
seems the most reasonable one. 

This choice can be infered using the scaling of the dressing functions extracted from 
the renormalization group equation, see ref. for details. The dressing functions can be 
expressed as 

where the running coupling provides the correct one loop scaling in the ultraviolet. Con- 
sequently the function R{x) approaches unity for high momenta. Furthermore, from the 
known infrared behavior of Z{x), G{x) and a{x) one infers that R{x) is proportional to x'^ 
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in the infrared. Writing Z\ in terms of a[x) and R[x) yields 



Zi{x,y,z;s,L) = 




In the perturbative region R{y),R{z) — > 1 and the function Zi is therefore slowly varying 
for any a and b according to the logarithmic behaviour of the running coupling a. In the 
infrared, however, a approaches its fixed point while the functions R behaves like a power. 
Consequently the choice a = b = 3S guarantees the weakest momentum dependence of Zi 
which is illustrated in Fig. Shown is the y-dependence of the function Zi{x, y,z = y; s, L). 
(Note also that due the symmetry Zi{x, y, z; s, L) = Zi{x, z, y; s, L) and the absence of an 
explicit x-dependence this is sufficient to demonstrate its momentum dependence.) In the 
perturbative momentum regime the function Zi does not vary with the parameter a. So 
all three choices give the same logarithmic running in momentum as required to give the 
correct one loop scaling behaviour of the integral. In the infrared, however, a change in a 
gives rise to substantial changes in the behaviour of Zi, with only the choice a = 35 leading 
to a constant. 



APPENDIX B: THE INFLUENCE OF ZERO MODES ON THE SOLUTIONS 

An important point when formulating the Dyson-Schwinger equations on the torus could 
be the treatment of the zero modes. In addition, on the torus an infrared analysis like the 
one in fiat Euclidean space-time is not possible, and one is left with the problem how the 
dressing functions behave at vanishing momenta. Guided by the intuition that especially 
the long ranged modes should be affected by the finite volume we assume in the following 
Z[x — 0) = just like in the continuum and G{x ^ 0) = const if zero modes are neglected. 
Phrased otherwise we assume that the zero modes are the missing ingredient to ensure the 
correct infinite volume limit for the torus results. Therefore, if on tori of different volumes 
G{x = 0) shows no sign of becoming divergent, the infrared enhancement seen in G{x — >■ 0) 
or in the flat space-time results has to be due to the torus zero modes of gluons and ghosts. 

Therefore, in this appendix, we will show that the assumption G{x = 0) < oo does not 
lead to a contradiction in the equations on the torus if zero modes are neglected. To this 
end we focus on the truncation scheme without angular approximations. First we rewrite 
equations (|T5|,|T6|) as 

9^N. I , (47) 



G{x) '' J (27r)4 xy 

1 ^ , 2^0 f d^q M{x,y,z) 

- Z-i + g — —— G[y)G[z) 



Z{x) 3 J (27r)4 xy 

f A Q(x,y,z)G(i/)-2-6^G(^)-2-65 



3 J (27r)4 xy Ziyf^ Z{z) 
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(48) 



According to Appendix A we have chosen a = b = 36, where S = —9/44, the anomalous 
dimension of the ghost. The kernels have the form: 

I f ix-yf\ 1 fx + y\ 1 sin^e 
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M{x, y,z) = - ——X + 77-7— + 7-' (50) 



2 4xy 2 2x 4x 

^ 1/1x3 19 -( 

Q[x,y,z) = + x — xy + ——y + - — 

z'^ \8 y 8 4 8x 

1 fx' 15 + C 17-C 
+- z — X — y + C— 

z \ y 4 4 X 

19-C x 3-4C 9Cy 
8 y 2 4 X 

+ (51) 

\x 4?/ / Sxy 

We first analyse the behavior of the integrands in the hmit ?/ — > for finite momenta x. 
Then Z{z) — > Z{x) and ^(2;) — G'(x) and the kernels times the respective dressing functions 
are to appropriate order in momentum y: 

mmK(,,,,,) ^ Gm(.f-^ , (52) 

xy x^ 

^M3^M{x,y,z) ^ G(0)G(x)l(l + (C-2)cos2e) , (53) 
xy x^ 

G{y)-'-''G{yy'''' _^ G(0)-^-^^G(x)-^-^^ / Ccos^9 \ 

xy \ xy J 

Furthermore, z = x ^ y — 2^/xycosG has been used and terms proportional to cos0 have 
been dropped, as they either integrate to zero in the continuum or cancel each other in the 
sums on the torus. Each of the expressions (p2| , |53| , p^ is then the appropriate term for j = 
on the right hand side of the Dyson-Schwinger equations on the torus. Clearly one observes, 
that only a finite ghost mode G(0) avoids trouble with divergencies. This is especially true 
for the kernel Q' of the gluon loop, as Z^^^{y 0) is more singular than the simple pole, so 
this kernel vanishes for small momenta y. The other two expressions (^) and ( p3|) are finite. 
One is then left with the ambiguous quantities sin^ G and cos^ which will be replaced by 
their integrals from zero to 27r in the calculation at the end of this section. The arbitrarines 
of this procedure is considerably milded by the observation that any number plugged in for 
the trigonometric functions yields the same qualitative result at the end of this section. 

Second, we take the limit z ^ 0, which on the torus is identical to G — > 0. The ghost 
kernel sin^ Q/ z' alone would certainly diverge as — >■ 0, but taking into account the power 
law behaviour Z{z) ~ z'^'^ for the gluon dressing function the integrand is zero in this 
limit. This is valid for k > 0.5, which is in agreement with the infrared analysis in the 
continuum. We therefore may omit the points 2; = in the ghost equation. The situation is 
different in the gluon equation where the kernel of the ghost loop has a finite limit 2; — > 0: 
M{x,x,0)/xy = (C + l)/(2x^). Therefore with a finite ghost dressing function G(0) the 
points z = in the ghost loop contribute but no divergencies occur. In the gluon loop the 
kernel Q' multiplied by the dressing functions approaches zero for vanishing momentum z 
due to the power law behaviour of Z~^^{z — » 0). 

To obtain a definite value for G{0) we now investigate the behaviour of the equations 
(^,^) in the limit x —>■ 0. The integrands are then given by 

^MMA-(x.y.^)-.Gfem,)2l4^. (55) 

xy y^ 
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FIG. 11: The ghostmode G{0) compared with the results for finite momentum x on the torus. 
For convenience we have kept the logarithmic momentum scale and plotted the zero modes on the 
left border of the figure. 



^ X, y, z) ^ G{y)G{y) ^ + ... , 56 

xy \ xy J 

where the kernel Q' is of no interest, as we know the gluon loop to be sub leading in the 
infrared. Clearly, the kernel M of the ghost loop in the gluon equation is now singular for 
X —>■ 0, corresponding to a vanishing gluon dressing function in the infrared. This result 
confirms our working hypothesis that the gluon mode Z{0) is not affected by the finite 
volume of the torus. The integrand of the ghost equation is finite up to the point y = 0. 
There the pole in the kernel is canceled by the behavior of the gluon dressing function 
Z{y) ~ y^'* resulting in a zero for vanishing momentum x and y. 

We therefore arrive at a consistent set of equations for a finite ghost mode G{0) and a 
vanishing gluon mode Z{0). In Fig. |TT| we show the results for the ghost dressing function 
gained on two different volumes on the torus. Within numerical accuracy the values of G{0) 
are the same for the two volumes. Obviously terms with high loop momentum y contribute 
most to the right hand side of the ghost equation for vanishing momentum x. Furthermore, 
one observes that the actual value of G{0) is not in accordance with an extrapolation of the 
ghost curves to the infrared. There is also no sizeable change of the gluon and the ghost 
dressing function when G{0) is set to zero by hand. This has been done in all calculations 
in the main body of this paper. 

Having shown that G{0) < oo on the torus even in the infinite- volume limit and assuming 
that the torus should provide a reasonable infrared regularization of physics in flat space- 
time we conclude that the divergence of G{0) is very probably due to the torus zero modes 
of gluons and ghosts. Noting furthermore that a diverging G{0) is related to the Zwanziger- 
Gribov horizon condition and the Kugo-Ojima confinement criterion this indicates a direct 
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relation between zero modes, the Gribov horizon and confinement. 



APPENDIX C: NUMERICAL METHODS FOR FLAT EUCLIDEAN SPACE- 
TIME 

The numerical methods to solve the coupled gluon-ghost system in flat Euclidean space- 
time after applying angular approximations has been described in great detail in refs. M, 



Note that different numerical techniques have been employed in refs. [0, and ref . . We 
have used all these techniques to verify the independence of our solutions from details of the 
numerical treatment. 

First, we discuss the numerical renormalization process. As we know the characteristics of 
the solutions a priori we may exploit the possibility of independent subtractions of the gluon 
and the ghost equat;ion. The procedure followed in ref. was to eliminate renormalization 
constants Z3 and Z3 by subtracting the equation at the renormalization scale s = The 
results presented in sect. IV. A are generated by solving the equations 



Z{x) Z[sz) 487r2 \ Jq x \ 2x J 2y 

Gfe) + 0G^(!/) j . (57) 



/o Sz V si. 2sz 
1 1 9g^N, ^ '■^ 



G[x) G(sg) 4 487r2 



Jo 3; X j^. y 



Z{sg) r ^^G{y) + r ^Z{y)G{y)\ , (58) 
io sgSg J so y 



for the truncation scheme of sect. 11. B and analogous equations for the scheme of sect. II. A 
thereby introducing subtraction points sz and sq- The parts of the integrals from y = to 
an infrared matching point e are carried out analytically, with the dressing functions in the 
integral beeing replaced by their leading infrared behaviour ( p2D 

Z{x) = Ax^", G{x) = Bx~\ 

The infrared analysis leads to flxed values k and ac = a(0) = {^^^)AB'^, so A is left as a 
free parameter in the infrared expansion. For numerical reasons it is convenient to choose 
the subtraction points sq = and to lie above the maximum of Z{x). The two input 
parameters that determine the renormalization point where we require the normalization 
conditions G{s = iJ?) = Z{s = /x^) = 1 are then A and the value the gluon dressing function 
at its subtraction point, i.e. Z{sz). 

The equations are solved using the Newton iteration method thereby generating values 
for the propagator functions. These are then plugged into the unsubtracted equations to 
obtain the respective values of the renormalization constants Z3 and Z3 for a given cutoff A 
and a given renormalization scale /i which allows flnally to determine the numerical solutions 
respecting the normalization conditions G{s = /i^) = Z{s = /x^) = 1. 
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